SPREADING OF LAGRANGIAN REGULARITY ON 
RATIONAL INVARIANT TORI 



JARED WUNSCH 

Abstract. Let Ph be a self-adjoint semiclassical pseudodifferential op- 
erator on a manifold M such that the bicharacteristic flow of the princi- 
pal symbol on T* M is completely integrable and the subprincipal sym- 
bol of Ph vanishes. Consider a semiclassical family of eigenfunctions, or, 
more generally, quasimodes Uh of Ph ■ We show that on a nondegenerate 
rational invariant torus, Lagrangian regularity of Uh (regularity under 
test operators characteristic on the torus) propagates both along bichar- 
acteristics, and also in an additional "diffractive" manner. In particular, 
in addition to propagating along null bicharacteristics, regularity fills in 
the interiors of small annular tubes of bicharacteristics. 



1. Introduction 

It is a well-known fact of semiclassical microlocal analysis, that the ana- 
logue of Hormander's theorem on propagation of singularities for operators 
of real principal type [7] holds for the semiclassical wavefront set (also known 
as "frequency set"): it propagates along null bicharacteristics of operators 
with real principal symbol [H [12]. Given a Lagrangian submanifold C of 
T*M, we may introduce a finer notion of regularity, the local Lagrangian 
regularity along C. We show here that on rational invariant tori in integrable 
systems, local Lagrangian regularity not only propagates along bicharacter- 
istics, but spreads in additional ways as well. 

Let Ph be a semiclassical pseudodifferential operator on a manifold M, 
with real principal symbol p (this is automatic if P is self-adjoint). Assume 
that the bicharacteristic flow of p is completely integrable. (In fact we only 
need to assume integrability locally, near one invariant torus.) Let Uh be a 
family of quasimodes of Ph, i.e. assume that \\{Ph — ^)uh\\L2 = 0{h^) for 
some G N, as /i J, either through a discrete sequence or continuously. 
(Note that this certainly includes the possibility of letting Uh be a sequence 
of actual eigenfunctions). Let C be an invariant torus in the characteristic 
set {p = A}. Then the bicharacteristic flow is by definition tangent to £, 
and we show (even in the absence of the integrability hypothesis) that La- 
grangian regularity propagates along bicharacteristics — this is Theorem A 
below. If a single trajectory is dense in £, then this is the whole story for 
propagation, as the set on which Lagrangian regularity holds is open, hence 
the whole torus either enjoys Lagrangian regularity or none of it does. At 
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the opposite extreme, if £ is a torus on which ah frequencis of the motion 
are rationahy related, we may ask the finer question: what subsets of the 
space of all orbits may carry Lagrangian regularity? The answer (assuming 
a nondegeneracy condition holds) turns out to be somewhat constrained: 
given a single orbit, Lagrangian regularity along a small tube around it im- 
plies Lagrangian regularity along the orbit itself. This is our Theorem B. 
(In the special case of two-dimensional tori, we can go further: again, either 
the whole torus enjoys Lagrangian regularity or no points on it do.) The 
order of regularity up to which our result holds is constrained by the order 
of the quasimode. We speculate that a finer theorem may be obtainable by 
more authentically "second-microlocal" methods. 

Example 1. As a simple example of our main result, Theorem B, we consider 
the case M = S^. x S^, Ph = h'^A = -h'^{d'^/dx^ + d^/dy^); we consider 
Lagrangian regularity on the Lagrangian torus £ = = 0, = 1} for 
quasimodes satisfying 

Lagrangian regularity on this particular C is special in that we may test for 
it using powers of the differential operator Dx = i~^{d/dx). The theorem 
tells us the following in this case: let T{x) be a smooth cutoff function 
supported on {|x| G [e,3e]} and nonzero at ±2e. Let cj) be another cutoff, 
nonzero at the origin and supported in [— 2e, 2e]. If, for all k' < k, we have 



{T{x)u,] 



< C < cx), 



then for all k' < k, 



< C < oo, 



D'J{^{x)uh) 

i.e. the regularity fills in the "hole" in the support of T. In this special 
case, the result can be proved directly by employing a positive commutator 
argument using only differential operators; the positive commutator will 
arise from the usual commutant h~^xDx- 

A less trivial example, that of the spherical pendulum, is discussed in ^ 
below. 

The methods of proof (and the idea of the paper) arose from work of Burq- 
Zworski [5, fS] and a subsequent refinement by Burq-Hassell-Wunsch [3] on 
the spreading of mass for quasimodes on the Bunimovich stadium. The 
central argument here is a generalization of the methods used to prove that a 
quasimode cannot concentrate too heavily in the interior of the rectangular 
part of the stadium (which is essentially the example discussed above on 
M = S^ X S^). 

We remark that our hypotheses in this paper are quite far from those 
in the study of "quantum integrable systems" where one examines eigen- 
functions of a system of n commuting operators on an n-manifold. For 
instance, if we take Ph = + h^V on the torus, with V a real valued. 
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smooth bump-function, then the operator satisfies the hypotheses of our 
Theorems A and B, and yet there does not exist a system of n — 1 other op- 
erators commuting with P^, with independent symbols. Moreover, even in 
the completely integrable case, given that we study eigenfunctions of a single 
operator, it may be possible to use the degeneracy of the system to construct 
non- or partially-Lagrangian quasimodes. Little seems to be known in this 
direction. 

The author is grateful to Andras Vasy for helpful discussions on La- 
grangian regularity, and to Clark Robinson for introducing him to isoen- 
ergetic nondegeneracy. He has also benefitted greatly from comments on an 
earlier version of the manuscript by Maciej Zworski, and by an anonymous 
referee. This work was supported in part by NSF grant DMS-0401323. 

2. LAGRANGIAN REGULARITY 

We begin by setting some notation and recalling some concepts of semi- 
classical analysis. For detailed background on this subject, we refer the 
reader to [Hill]. 

Let be a smooth manifold and fix £ C T*M a Lagrangian subman- 
ifold. Throughout the rest of the paper, we assum^l E L'^[M\Q}/'^), 
with h G (0, /lo); here Vt^/"^ denotes the bundle of half- densities on M, 
i.e. the square root of the density bundle lA^Mj. We will in future, how- 
ever, suppress the half-density nature of Uh as well as its /i-dependence, 
writing simply u € L'^{M); similarly, all operators will tacitly be semi- 
classical families of operators, operating on half-densities. The hypothesis 
that our operators act on half-densities ensures that if ^ = Op/j(a) with 
a(x,^; h) ~ ao(x,^) -|- /iai(x,^) -|- . . . , the terms ao (principal symbol) and 
ai (subprincipal symbol) are both invariantly defined as functions on T*M 
(see m)- 

Furthermore, we will deal with an operator P rather than P—X, absorbing 
the constant term into the definition of the operator. 

We begin by defining a notion of Lagrangian regularity of a family of 
functions along £, following the treatment of the "homogeneous" case in [dj- 

Definition 2. Let A4 denote the module (over ^'/^(M)) of semiclassical pseu- 
dodifferential operators with symbols vanishing on C 

Let q £ C, k £ N, and u £ L?'{M). We say that u has Lagrangian regularity 
of order k at q, and write q G S^{u), if and only if there is a neighborhood 
[/ of g in T*M such that for all fc' = 0, 1, . . . , and all Ai, . . . , At' £ M 
with WF' Aj C U, h-^' Ai ■ ■ ■ Ak>u £ L^{M). 

Proposition 3. Fix q £ C, and let Ai (i = l,...,n) be a collection of 
elements of Ai with da{Ai) spanning N*C. We have 

q £ Sc{u) h^'"' Ai^- ■ ■ Ai^,u £ L'^ V(zi, . . . , ifc/) e {1, . . . , n}''', k' < k. 



We may just as well assume that h I through a discrete sequence; this will make no 
difference in what follows. 
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Proof. We begin with the case k = 1. Given any B characteristic on C and 
microsupported sufficiently close to q, we may factor a{B) = Y^Cia{Ai) by 
Taylor's theorem. Thus, letting Cj be operators with symbol q, we obtain 

h-^Bu = h-^dAiU + Ru 

for some semiclassical operator hence we obtain the desired estimate on 
h~^Bu since R is uniformly (in h) L^-bounded. 

More generally, if Ba^ , • • • , -Bj^. is a A;-tuple of operators characteristic on 
£, we have 

k 

h-^Bi, ■ ■ ■ Bi^u = \[{Ci^Ai. + hRi^)u; 
i=i 

We then obtain the desired estimate inductively, using the fact that each 
commutator of the form [C, A\ or [i?, A\ produces a further factor of h. □ 

We note that it follows from the work of Alexandrova [Ij that S'^ = £ if 
and only if we can actually write u in the form of an oscillatory integral 

plus a term with semiclassical wavefront set away from £; here i;^ is a 
phase function parametrizing the Lagrangian C in the sense introduced 
by Hormander. This is the semiclassical analog of a central result in the 
Hormander-Melrose theory of conic Lagrangian distributions [9, Chapter 
25]. 

We now observe that the analogue of Hormander's theorem on propaga- 
tion of singularities for operators of real principal type is easy to prove in 
our setting. 

Theorem A. Let P G ^'/^(M) have real principal symbol p. Let £ C {p = 0} 
he a Lagrangian submanifold ofT*M. Then Pu G h^~^^L'^{M) implies that 
S^{u) is invariant under the Hamilton flow of p. 

The author is grateful to M. Zworski for suggesting the following brief 
proof. 

Sketch. By [HI Theorem 21.1.6], there is a local symplectomorphism taking 
p to ^1 and £ to £o = {'^ = 0}- Following the development in p], we 
may quantize this to a semiclassical FIO that conjugates P to hD^i modulo 
0{h°^) (cf. [9, Theorem 26.1.3] in the non-semiclassical setting) . Lagrangian 
regularity along Cq is iterated regularity under h~^{hD^i), i.e. is just classical 
Sobolev regularity, uniform in h. The theorem thus reduces to the statement 
that Sobolev regularity for solutions to D^iu G h^L'^{M) propagates along 
the lines (x^ G M, = const), which is easily verified. □ 
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3. Integrable flow 

We continue to assume that P G ^/^(M) has real principal symbol. We 
now further assume that p = cr(P) has completely integrable bicharacteristic 
flow, i.e. that there exist functions /2, . . . , /n on T*M, Poisson commuting 
with p and with each other, and with dp, df2, • • • , dfn pointwise linearly in- 
dependent. We again emphasize that we in fact only require the /j's to 
exist in some open subset of interest in T*M. Let S denote the characteris- 
tic set in T*M. Let {Ii, . . . , In,Oi, . . . ,9n) be action-angle variables and let 
LJi = dp/dii be the frequencies. We also let Uij = d^I/dlidlj. (We refer 
the reader to [2] for an account of the theory of integrable systems, and in 
particular for a treatment of action-angle variables.) 

Let £ C S be a rational invariant torus, i.e. one on which Ui/iVj £ Q 
for all i,j = We further assume that C is nondegenerate in the 

following sense: we assume that the matrix 

fuJn ... iVin UJl\ 

(1) 

UJni . . . LOnn 

\^ ui ... y 

is invertible on C This is precisely the condition of isoenergetic nondegen- 
eracy often used in KAM theory (see [2], Appendix 8D). It is easy to verify 
that the condition is equivalent to the condition that the map from the 
energy surface to the projectivization of the frequencies 

{p = 0} 9 / ^ [uji{l) ■ ■ ■ : ujn{I)] G MP" 

be a local diffeomorphism. 

For later convenience, we introduce special notation for the frequencies 
and their derivatives on £ : we let 

On C, we of course only know from Theorem A that S'^{u) is a union of 
orbits of Hp, which, being rational, are not dense in C. There are, however, 
further constraints on S'^{u). 

Definition 4. An annular neighborhood of a closed orbit p is an open set 
U = V\K C C such that p C K C V with K compact and V open in C 

We can now state our main result. 

Theorem B. Suppose Pu = / G h^'^^L?'. Let p be a null bicharacteristic for 
p on the rational invariant torus C If a small enough annular neighborhood 
of p is in S^{u), then so is p. 

The meaning of "small enough" depends only on the uJi 's. 

If n = 2 then either Sc{u) = C or Sc{u) = 0. 

Thus, conormal regularity propagates "diffusively" to fill in annular neigh- 
borhoods. 
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Example 5. Horozov [10^ [TT] has studied the spherical pendulum, i.e. the 
system on T*S^ with Hamiltonian h = {l/2)\^\ +x-^ onT* (with one of 
the Euchdean coordinates on 5^^ C M^). Integrals of motion are h and pg, the 
angular momentum. Horozov showed that when h E (—1,1] U [7/\/T7, oo), 
all values of pg lead to isoenergetically nondegenerate invariant tori, while 
for /i G (1, 7/\/T7), there are exactly two values of pe for which isoenergetic 
nondegeneracy fails. Thus our results show that if we consider quasimodes 
for the operator 

Ph = {l/2)h^As2 + X3 

then for any torus £ not associated to the one of the exceptional pairs of 
{h,pg) identified by Horozov, either Sc = CoySc = 0- 

Example 6. We now illustrate with an example the necessity of the isoener- 
getic nondegeneracy condition. As in the introduction, let M = x , but 
now let P = hDx] it is easy to verify that no Lagrangian torus is isoenerget- 
ically nondegenerate in this case. Let £ = = 77 = 0}, the zero-section of 
T*M. Lagrangian regularity in this setting is, as noted above, just Sobolev 
regularity, uniform in h. 

Let "0(2/) be a bump function supported near y = 0. Then 

u{x,y) = e^^(^)/^ 

has wavefront set only in £. It is manifestly Lagrangian on the complement 
of suppV', which forms an annular neighborhood of the orbit {x ^ S^,y = 
0, ^ = = 0} C supp^/'. It is not Lagrangian, however, on suppV', as it lacks 
iterated regularity under h~^{hDy). 

4. Symbol Construction 

By shifting coordinates, we may assume that p is the orbit passing through 
{9 = 0}. 

For each i,j let 

jij = minkj(zz{{Oi + 27rk)ujj - {Oj + 2TTl)uJi), 

where min denotes the value with the smallest norm, i.e. may be positive 
or negative. Each jij then takes values in an interval determined by tJj, 
iJj, and is smooth where it takes on values in the interior of the inter- 
val. (If uJi = p/q and Uj = p' /q' then ^ij takes values in [— 7ra,7ra] where 
a = gcd{qp' ,pq')/qq' .) The "small enough" condition in the statement of 
Theorem B is just the following: each should be smooth on the annular 
neighborhood of p where we assume Lagrangian regularity. 

Note that ^ij{0) = for all i,j exactly when there exists G R", equiv- 
alent to 9 modulo 27rZ" , such that [61 : . . . On] = \uJi : ■ ■ ■ : Thus the 
functions define p on £ : we have {/ = /, 7,^- = Vi,j} = p; indeed, 
the vanishing of each 7i,i+i and of 7^,1 suffices to define p, and these n 
functions may be taken as coordinates on £ in a neighborhood of p. The 
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central point of our argument will be that the jij are "propagating vari- 
ables" with derivatives along the flow that, taken together, will suffice to 
give Lagrangian regularity. 

Since the jki define p and are smooth on the annular neighborhood U 
where we have assumed regularity, there is a smooth cutoff function 

V' := ^"(712,, 723, • • • ,7n-l,n,7n,l) 

with = 1 on p and Vip having its support on C contained in U. We may 
also arrange for ip to be the square of a smooth function. Let (j)^: be a cutoff 
supported in [— e,e], with smooth square root. 
Let 

(2) aij{x) = ■ ■ ^ij ■ {iViUJj - OJjUJi) 

We compute first that, where 'yij G C°°, 

(3) {p, 7ij} = {uJiUj - Lopi) 

(since is locally given by expressions of the form [{Oi + 2T:k)ujj — {9j + 
2Trl)iJi) with k,l fixed) and hence that 

{p,aij} = {p,tp}<j)e{\l - ^IhijiuJiuJj -uJjUi) + ip(j)e{\l - l\){ujiujj -ujjUif. 

We further note that as V' is a function of the jij^s, by ([3|) the first term in 
this expression is a sum of terms divisible by 

for various ki,li. Thus we may write 

(4) {p, ttij} = ^ Ckfl + ^^(peill - 7\){uJiUJj - UJjUJif, 

where each and fi vanishes on C and with support intersecting C only in 
U. 

We will also employ a symbol that is invariant under the flow: for each 
j = 1, . . . ,n, set 

Wj = - 
5. NONDEGENERACY 

Using a positive commutator argument, we will flnd that we can control 
operators whose symbols are multiples of (wjcUj — ujjUJi). These quantities 
vanish on i2, but our nondegeneracy hypothesis permits us to use them to 
control Lagrangian regularity on C To see this, rewrite 

{ijOiUJj — uJjoJi) = {uJi — oJi)Zoj — {ojj — TjOj)uJi 

and expand about C in the / variables, to rewrite this as 

y^(uJifcuJj - TDjkU}i){Ik - Ik) + 0{{I - 7)^) 

k 

We now prove a key algebraic lemma: 
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Lemma 7. Let vi, . . . ,Vn, and Vij, i,j = l,...,n be real numbers, with 
Vij = Vji. The functionals aij{x) = Ylki'^ikVj - VjkVi)xk (for i,j = 1,... ,n) 
together with the covector (vi, . . . ,Vn) span {R"')* iff the matrix 



(5) 



is nondegenerate. 



( v\\ 



'-'nn '-'n 
V„. 



Proof of Lemma. We may assume that not all of the ViS are zero, as the 
result is trivial in that case. 
Let 

/ VliVj - VljVi 

Cij = : 
Letting A be the matrix with entries Vij and 



where e,- is the standard basis for 



Cr 



we have 



Auij . 



Let U denote the span of the Ujj's. Thus, 



where v = {vi, . . . Vn) ■ Thus, U = v . Hence the span of the Qj is of A{v ). 
The assertion of the lemma is then that A{v-^) and v are complementary iff 
the matrix ([5]) is nondegenerate. This equivalence follows from the observa- 
tion that 



A V 




AiS + zv 
{v,w), 



hence ([5]) has nontrivial nullspace iff there exists a nonzero w ^ v-^ with 
Aw G Rv. □ 



6. Proof of Theorem B 

We note, first of all, that in the special case when n = 2, a neighborhood 
of any closed orbit p' ^ p is itself an annular neighborhood of p. Hence the 
special result for n = 2 follows directly from the general one. 

We now prove Theorem B by induction on k; we suppose it true for 
k < K — 1 (and note that for = it is vacuous). 
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Let Aij G ^'/i(M) be self-adjoint, with symbol Uij constructed above and 
vanishing subprincipal symbol. Then we have by ([4]), 

(6) ih-^[P, A,] = h-'^Bl + h-'^EkFi + R 

k,l 

with Bij self-adjoint with vanishing subprincipal symbols, and 

(7) (y{Bij) = bij = (^/;</>e(|/ -7|))^/^ • {uJiUJj -UJjUJi), 

and with Ef^, Fi characteristic on C with the supports of alE/.), cr(Fi) inter- 
secting C only in U. {R, Ek, and Fi of course depend on i,j but we suppress 
these extra indices.) 

Let Wj have symbol Wj constructed above, and be self-adjoint with van- 
ishing subprincipal symbol. Then 

ih-^[P,Wj] G ^h{M). 

For a multi- index a with \a\ = K — 1, set 

We will also need the operator denoted in multi-index notation 

W"" = W^' . . . Wn''-' 

Now we examine 

(8) ih-^''-^{{P*Qi, - Qi,Pu,u) = ih-^^'-'mjuJ) - {f,Q*ju)). 
For any 5 > 0, we may estimate the RHS by 

C5\\h-^-'ff + S{\\h~^<Q,juf + \\h-^Q*^uf). 

Note that both Qij and Q*j are (2K+l)-fold products of operators vanishing 
on M, and that each contains the factors Aij and W". By ^ and ([7]), a{Aij) 
is divisible by a{Bij); thus, by elliptic regularity we may estimate the RHS 
by 

2K-2 

Cs\\h'^~^ff + C6j2\\f^'^BijW''uf + h~^{DjU,u) 

a j=0 

where C is independent of S, and each Dj is a sum of products of j elements 
of Ai, all microsupported on suppAj^-; these arise from commutator terms 
in which we have reordered products of elements of Ai. 

Now we recall that P* - P e h?-^h{M), hence, by ^ we may write 

(9) {ih-^^-\P*Qi, - QijP)u,u) 

2K-2 

= \\h-^'^BijW''uf + Y^{h-^EkW''u,h~^FiW''u)+ ^ h~^(DjU,u^ 

k,l j=Q 

with the Dj sharing the properties of the Dj above. 
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Putting together the information from our commutator, we now have, for 
all 5>0, 

(10) {l-C6)\\h~^^BijW"u\f 

2K-2 

< Cs\\h-^-^ff+Y.^\\h~^EkW''uf+\\h-^FiW''uf)+ ^ h'^ 



DjU, u 



with the -Dj's satisfying the same properties as Dj above. Each of the 

and Fi terms is controlled by our hypothesis on S^u, while the Dj terms are 
bounded by the inductive assumption. 

Now we use our nondegeneracy hypothesis as reflected in Lemma[7l Recall 
that £ C E, hence the operator P is characteristic on £; moreover, we have 
dp\c = Yl^kdlk, hence Lemma [7] tells us that P and Bij, for i, j = 1, . . . ,n, 
are a collection of operators fulfilling the hypotheses of Proposition [3l Thus, 
adding together equations (fTUj) for all possible values of and multi-index 
a, together with terms involving P rather than Bij (which vanish up to 
commutators of P with VF's), we obtain the desired estimate, by Proposi- 
tion El □ 
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